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Abstract 

We investigate the particle trajectories in a constant vorticity shallow water flow over a 
flat bed as periodic waves propagate on the water's free surface. Within the framework 
of small amplitude waves, we find the solutions of the nonlinear differential equations 
system which describes the particle motion in the considered case, and we describe the 
possible particle trajectories. Depending on the relation between the initial data and 
the constant vorticity, some particle trajectories are undulating curves to the right, 
or to the left, others are loops with forward drift, or with backward drift, others can 
follow some peculiar shapes. 

1 Introduction 

The motion of water particles under regular waves which propagate on the wa- 
ter's free surface is a very old problem. It was widely believed that the particle 
trajectories are closed. After the linearization of the governing equations for 
water waves, analysing the first approximation of the nonlinear ordinary dif- 
ferential equations system which describes the particle motion, one obtained 
that all water particles trace closed, circular or elliptic, orbits (see, for exam- 
ple, Q3], [3D], [22], [25, [23, [23, [2H]-a conclusion apparently supported by 
photographs with long exposure [14], [27], [28]). 

While in this first approximation all particle paths appear to be closed, 
in [llj it is shown, using phase-plane considerations for the nonlinear system 
describing the particle motion, that in linear periodic gravity water waves no 
particles trajectory is actually closed, unless the free surface is flat. Each par- 
ticle trajectory involves over a period a backward/forward movement, and the 
path is an elliptical arc with a forward drift; on the flat bed the particle path 
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degenerates to a backward/forward motion. Similar results hold for the parti- 
cle trajectories in deep-water, that is, the trajectories are not closed existing a 
forward drift over a period, which decreases with greater depth (see |5J). These 
conclusions are in agreement with Stokes' observation [33]: "There is one result 
of a second approximation which may possible importance. It appears that the 
forward motion of the particles is not altogether compensated by their backward 
motion; so that, in addition to their motion of oscillation, the particles have a 
progressive motion in the direction of the propagation of the waves. In the case 
in which the depth of the fluid is very great, this progressive motion decreases 
rapidly as the depth of the particle considered increases." 

For linearized irrotational shallow water waves, there are obtained very recently 
in [19] the exact solutions of the nonlinear differential equations system which 
describes the particle motion. Beside the phase-plane analysis, the exact solu- 
tions allow a better understanding of the dynamics. In [115] it is shown that 
depending on the strength of underlying uniform current, beneath the irrota- 
tional shallow water waves some particle trajectories are undulating path to the 
right or to the left, some are looping curves with a drift to the right. 
The steady (traveling) linear water waves with constant vorticity were studied 
recently in |15j . |16j . Here linearity means that the waves are small perturba- 
tions of shear flows. The linear system obtained in this way is solvable. Further, 
making a phase portrait study in steady variables for the nonlinear differential 
equations system which describes the particle paths, it is found that for positive 
vorticity, the steady wave resembles that of the irrotational situation, though for 
large enough vorticity the particles trace closed orbits within the fluid domain. 
For negative vorticity all the fluid particles display a forward drift. 
Using the same approach as in |19j within the framework of linear water waves 
theory, we investigate in this paper the particle trajectories beneath shallow 
water waves in a constant vorticity flow. The obtained results and the way in 
which the paper is organized are presented below. 

Let us now give some references on the results obtained for the govern- 
ing equations without linearization. Analyzing a free boundary problem for 
harmonic functions in a planar domain, in [5] it is shown that there are no 
closed orbits for Stokes waves of small or large amplitude propagating at the 
surface of water over a flat bed; for an extension of the investigation in [5] to 
deep-water Stokes waves see [18]. Within a period each particle experiences a 
backward/forward motion with a slight forward drift. In a very recent preprint 
[TO] , the results in [5] are recovered by a simpler approach and there are also 
described all possible particle trajectories beneath a Stokes wave. The particle 
trajectories change considerably according to whether the Stokes waves enter a 
still region of water or whether they interact with a favorable or adverse uni- 
form current. Some particle trajectories are closed orbits, some are undulating 
paths and most are looping orbits that drift either to the right or to the left, 
depending on the underlying current. 

Analyzing a free boundary problem for harmonic functions in an infinite planar 
domain, in [7] it is shown that under a solitary wave, each particle is trans- 
ported in the wave direction but slower than the wave speed. As the solitary 
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wave propagates, all particles located ahead of the wave crest are lifted while 
those behind have a downward motion. 

Notice that there are only a few explicit solutions to the nonlinear governing 
equations: Gerstner's wave (see [T7] and the discussion in [3])), the edge wave 
solution related to it (see [1]), and the capillary waves in water of infinite or 
finite depth (see [13], [23] )■ These solutions are peculiar and their special fea- 
tures (a specific vorticity for Gerstner's wave and its edge wave correspondent, 
and complete neglect of gravity in the capillary case) are not deemed relevant 
to sea waves. 

The present paper is organized as follows. In Section 2 we recall the govern- 
ing equations for gravity water waves. In Section 3 we present their nondimen- 
sionalisation and scaling. It is natural to start the investigation for shallow water 
waves by simplifying the governing equations via linearization. The linearized 
problems for an irrotational shallow water flow and for a constant vorticity shal- 
low water flow are written in Section 4. We also obtain here the general solutions 
of these two linear problems. In the next section we find the solutions of the 
nonlinear differential equations systems which describe the particle motion in 
the two cases, and we describe the possible particle trajectories beneath shallow 
water waves. We see that these particle trajectories are not closed. Section 
5.1 contains the irrotational case. Depending on the strength of the underlying 
uniform current, the particle trajectories are undulating path to the right or to 
left, are looping curves with a drift to the right, and, if there is no underlying 
current or the underlying current is moving in the same direction as the irrota- 
tional shallow water wave with the strength of the current smaller than 2, then, 
the particle trajectories obtained are not physically acceptable (Theorem 5.1). 
In dealing with the linearized problem not with the full governing equations, 
we expect to appear solutions which are not physically acceptable. Section 5.2 
contains the case of a constant vorticity flow. Also in this case the particle 
trajectories are not closed. Depending on the relation between the initial data 
(xo, z Q ) and the constant vorticity cjq, some particle trajectories are undulating 
curves to the right, or to the left, others are loops with forward drift, or with 
backward drift, others can follow some peculiar shapes (Theorem 5.2, Theorem 
5.3). 



2 The governing equations for gravity water waves 

We consider a two-dimensional inviscid incompressible fluid in a constant grav- 
itational field. For gravity water waves these are physically reasonable assump- 
tions (see [20] and [25]). Thus, the motion of water is given by Euler's equations 

u t + uu x + vu z = ~- p Vx , 
v t + uv x + vv z = -- p Pz - g 

Here [x, z) are the space coordinates, (u(x, z, t), v(x, z, t)) is the velocity field of 
the water, p(x, z, t) denotes the pressure, g is the constant gravitational acceler- 
ation in the negative z direction and p is the constant density. The assumption 
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of incompressibility implies the equation of mass conservation 

u x + v z = (2) 

Let ho > be the undisturbed depth of the fluid and let z = ho+rj(x, t) represent 
the free upper surface of the fluid (see Figure 1). The boundary conditions at 
the free surface are constant pressure 

p = pa on z = h + r)(x, t), (3) 

Po being the constant atmospheric pressure, and the continuity of fluid velocity 
and surface velocity 

v = Tj t + ur\ x on z = h Q + rj(x, t) (4) 

On the flat bottom z — 0, only one condition is required for an inviscid fluid, 
that is, 

d = oaz = (5) 

Summing up, the exact solution for the water-wave problem is given by the 
system (JTJ)— (JSJ) - In respect of the well-posedness for the initial- value problem 
for (QJ-lIS]) there has been significant recent progress, see [12] and the references 
therein. 

A key quantity in fluid dynamics is the curl of the velocity field, called 
vorticity. For two-dimensional flows with the velocity field (u(x, z, t), v(x, z, t)), 
we denote the scalar vorticity of the flow by 

u(x,z) = u z - v x (6) 

Vorticity is adequate for the specification of a flow: a flow which is uniform 
with depth is described by a zero vorticity (irrotational case) , constant non-zero 
vorticity corresponds to a linear shear flow and non-constant vorticity indicates 
highly sheared flows. See in the Figure 1 an example of a linear shear flow with 
constant vorticity u> = const := luo > 0. 
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Figure 1. Linear shear flow. 



The full Euler equations CO)-© are often too complicated to analyze directly. 
One can pursue for example a mathematical study of their periodic steady 
solutions in the irrotational case (see see [T], [30 ) or a study of their periodic 
steady solutions in the case of non-zero vorticity (see [8], [9]). But in order to 
reach detailed information about qualitative features of water waves, it is useful 
to derive approximate models which are more amenable to an in-depth analysis. 
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3 Nondimensionalisation and scaling 

In order to develop a systematic approximation procedure, we need to char- 
acterize the water-wave problem ([I])-© in terms of the sizes of various fun- 
damental parameters. These parameters are introduced by defining a set of 
non-dimensional variables. 

First we introduce the appropriate length scales: the undisturbed depth of wa- 
ter ho, as the vertical scale and a typical wavelength A (see Figure 1), as the 
horizontal scale. In order to define a time scale we require a suitable velocity 
scale. An appropriate choice for the scale of the horizontal component of the 
velocity is y/gho- Then, the corresponding time scale is ^ h and the scale for 

the vertical component of the velocity is h , The surface wave itself leads 
to the introduction of a typical amplitude of the wave a (see Figure 1). For 
more details see [3D]. Thus, we define the set of non-dimensional variables 

x i — ^ Xx, z M> hnz, r\ H- an, t i-» , A , t, 

u i y y'ghou, v i-> ho v 

where, to avoid new notations, we have used the same symbols for the non- 
dimensional variables x, z, r), t, u, v, on the right-hand side. The partial 
derivatives will be replaced by 



ghp „, „, ^ \Zghg „, qi ^ ^gh „ 

(8) 



Vt 2$Vt, V X M- h ^V x , V Z M- &K Vz 



Let us now define the non-dimensional pressure. If the water would be station- 
ary, that is, u = v = 0, from the equations ^ and ([3]) with 77 = 0, we get for 
a non-dimensionalised z, the hydrostatic pressure p + pgh (l — z). Thus, the 
non-dimensional pressure is defined by 

P >-> Po + pgho(l - z) + pgh p (9) 

therefore 

Px i-> P^-Px, Pz >-> -pg + pgPz (10) 

Taking into account ©, ©, © and (fTU]), the water-wave problem QJ-© 
writes in non-dimensional variables, as 

u t + uu x + vu z = -p x 
5 2 (v t + uv x + vv z ) = -Pz 

Ux+V z =0 (11) 

v = e(r/ t + urjx) and p = en on z = 1 + erj(x, t) 
v = on z — 

where we have introduced the amplitude parameter e = and the shallowness 
parameter 5 = In view of ([Sjl, the vorticity equation (J5J writes in non- 
dimensional variables as 

r2„. , Vi^O 



= 8 A v x + z) (12) 
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For zero vorticity flows (irrotational flows) this equation writes as 

u z = 6 2 v x (13) 

For constant non-zero vorticity flows, that is, U)(x, z) =const:=wo, the equation 
(fH?]) becomes 

u z = o v x -\ lu q (14) 

9 

After the nondimensionalisation of the system ((I])-© let us now proceed 
with the scaling transformation. First we observe that, on z = 1 + erj, both 
v and p are proportional to e. This is consistent with the fact that as e — > 
we must have v — » and p — > 0, and it leads to the following scaling of the 
non-dimensional variables 

p i->- ep, (u, v) h-t e(it, v) (15) 

where we avoided again the introduction of a new notation. The problem (jllj) 
becomes 

u t + e(uu x + vu z ) = -p x 
6 2 [v t + e(uv x + vv z )] = -p z 

u x + v z = (16) 
v = rjt + £ur] x and p = r\ on z = 1 + erj{x, i) 
v = on z — 



(17) 



and the equation (fT2|) keeps the same form. 

In what follows we will consider in turn the cases of an irrotational flow and a 
constant vorticity flow. The system which describes our problem in the irrota- 
tional case is given by 

Ut + e(uu x + vu z ) = -p x 
8 2 [v t + e(uv x + vv z )} = -p z 
u x + v z = 
u z = 6 2 v x 

v = rjt + £ur] x and p — r\ on z = 1 + (Tj(x, i) 
v = on z = 

In the constant vorticity case, the problem is described by the following system 

u t + e(uu x + vu z ) = -p x 
S 2 [v t + e(uv x + vv z )] = -p z 
u x + v z = 
u z = S 2 v x + 



2„, i VSS,.,„ ( 18 ) 
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v = rjt + eu^x and p — i] on z = 1 + £77(2:, f) 
v = onz = 
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4 The linearized problem 

The two important parameters e and S that arise in water-waves theories, are 
used to define various approximations of the governing equations and the bound- 
ary conditions. The scaled versions (fTT)) and (fT5)l of the equations for our prob- 
lem, allow immediately the identification of the linearized problem, by letting 
e — > 0, for arbitrary S. The linearized problem in the shallow water regime is 
obtain by letting further 8 — > 0. Thus, we get the following linear systems, in 
the irrotational case 

u t + p x = 

Pz =0 
u x + v z = 

u z = 

v = rj t and p = r\ on z = 1 
u = on z = 

and in the constant vorticity case 

Ut + Px = 

+ w z = 

w z = ^7-^0 
v = r}t and p ^ rj on z = 1 
u = on z = 



(19) 



(20) 



From the second equation in (|19|). respectively (|20|) . we get in the both cases 
that p does not depend on z. Because p = rj(x, t) on z — 1, we have 

p = rj(x,t) for any < z < 1 (21) 

Therefore, using the first equation and the fourth equation in (|19p . respectively 
(|2T)|) . we obtain, in the irrotational case 



= I Vx(x, s)ds + T(x) 
Jo 



(22) 



and in the constant vorticity case 



u = — f i] x (x, s)ds + F(x) H — ° z (23) 
Jo 5 

where T is an arbitrary function such that 

T{x) = u(x, 0, 0) (24) 



Differentiating (|2^|) . respectively (1231) . with respect to a; and using the third 
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equation in (fT!?)) . respectively (|2T)j) . we get, after an integration against z, the 
same equation in the both cases 



v = —zu x = 



f Vxx{x, s)ds - J"' (a;) J 



(25) 



In view of the fifth equation in ([19]) , (|20[) . we get after a differentiation with 
respect to t, that r\ has to satisfy the equation 

Vtt - Vxx = (26) 

The general solution of this equation is i](x,t) — f(x — t) + g(x + t), where / 
and g are differentiable functions. It is convenient first to restrict ourselves to 
waves which propagate in only one direction, thus, we choose 

r,(x,t)=f(x-t) (27) 

From ([25]) . (|27|) and the condition v = r\t on z = 1, we obtain 

T(x) = f{x) + c (28) 

where Co is constant. 

Therefore, in the irrotational case, taking into account (|2"TT) . ([HI), (|25p . (|2T|) and 
(|28|) . the solution of the linear system (Tl9t is given by 

f?0M) = f{x-t) 

p(x,t) = f(x-t) , 2Q , 
u(x, z, t) = f(x - t) + c 

z, t) = —zf'{x — t) = —zu x 

In the case of a constant vorticity flow, from ([2"I]). (|2^|) . (gSJ, (HZJI and pg]). the 
solution of the linear system ([20)) is given by 

V( x ,t) = f(x-t) 
p(x, t) = f(x — t) 

u(x,z,t) = f(x-t) + ^^z + c (30) 
v(x, z, t) = —zf'(x — t) = —zu x 

5 Particle trajectories beneath linearized shallow water waves 

Let (x(t), z(t)) be the path of a particle in the fluid domain, with location 
(x(0), z(0)) at time t = 0. The motion of the particle is described by the 
differential system 

1 = ^,2,4) ldlJ 
with the initial data (x(0), z(0)) := (a^Oj^o)- 
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5.1 The case of an irrotational flow 

This case was investigated in [TH] . We summarize below the results presented in 
detail in reference [19] . Some new observations concerning the solutions shown 
in Fig. 2 and Fig. 6 are added. 
Making the Ansatz 

f{x -t) = cos(2tt(x - t)) (32) 
from (|2"9"|) . the differential system (|3"Tj) becomes 

|| = cos(2tt(:e - t)) + c 

(33) 

*g = 27rzsin(27r(a; - t)) 

Notice that the constant cq is the average of the horizontal fluid velocity over 
any horizontal segment of length 1, that is, 



1 



a:+l 



Co = — / u(s,z,t)ds, (34) 

1 Jx 

representing therefore the strength of the underlying uniform current. Thus, 
c = will correspond to a region of still water with no underlying current, 
Co > will characterize a favorable uniform current and Co < will characterize 
an adverse uniform current. 

The right-hand side of the differential system (|33|) is smooth and bounded, 
therefore, the unique solution of the Cauchy problem with initial data (xo,zq) 
is defined globally in time. 

To study the exact solution of the system (|3"3"|) it is more convenient to re- 
write it in the following moving frame 

X = 2ir{x-t), Z = z (35) 

This transformation yields 

^ = 27rcos(X) +2tt(co - 1) 

(36) 

H = 27rZsin(X) 
For Co = 0, we get the following exact solution of the system (|36l) : 

X(t) = 2arccot (27rf + a), a — constant (37) 

-f + {2irt + a) 2 - S 



Z(t) = Z(0) exp ( J 27rsin(X(s)) ds j = Z(0) exp ( In 



1 



(38) 
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Taking into account (|55|) . (|57|) and (1551) . we obtain the solution of the system 
(f33l) with the initial data (a;n, 2n): 

x(t) =t+ iarccot (2nt + a) 

(39) 

z (t) = T f^[l + (2irt + a) 2 ] 
From the initial conditions, we have a := cot (ttxq)- 

Studying the derivatives of x(t), z(t) with respect to t and the limits of x(t), 
z(t), for t — > ±oo, we draw the graph of the parametric curve (|39[) : 



Figure 2. Particle trajectory, in the irrotational case, for Co=0 



As is shown in Figure 2, the coordinate z increases indefinitely in time, thus, 
the solution (|39| is not physically acceptable. In dealing with the linearized 
problem not with the full governing equations, we expect to appear solutions 
which are not physically acceptable. 



For Co(co — 2) > 0, the exact solution of the system (|5fj|) has the following 
expression: 



X(t) = 2arccot <£n tan (a(t)) 



(40) 



Z(t) = Z(0)exp (j - 



47r£o tan (a(s)) 



rds 



Co tan (a(s))) 



(41) 



where 



CO 



a (t) ~ ~^±(2Trt + a) 



(42) 



(43) 



From 



(|4H)) and (HU, we obtain the solution of the system with the 
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initial data (xo,Zq), zq > 0, 

x{t) = t+ iarccot £q tan (a(f )) 



Z(t) = Z exp /„' 4 ^0tan(a(,)) , 

\ 1+ Co tan(a(s))) / 



(44) 



Studying the derivatives of x(t) and z(t) with respect to f, we get the following 
graphs: 



Figure 3. Particle trajectory, in the irrotational case, for c < - 1 




Figure 4. Particle trajectory, in the irrotational case, for - 1< c < 



Figure 5. Particle trajectory, in the irrotational case, for C(f> 2 

For Co G (0, 2], the exact solution of the system (|36|) has for | cot (=§-) | > Mo the 
following expression: 



X(f)=2arccot Mo coth (/3(f)) , Z(t) = Z(0)exp 

and for | cot (^) | < Mo the following expression: 
X(f) = 2arccot tanh (/3(f)) , Z(f) = Z(0)exp 



AttRq coth(/3(s)) 



T ds , 



J?o coth (/3(s))) 
(45) 



47rJ?otanh(/3(s)) , N 

— : 2^ s 

Mo tanh (/3(s))) / 

(46) 
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In the above formulas we have denoted by 



£n := 



2 -c 



Co 



m := £°*° (27 rt + a) 



(47) 



(48) 



From (|55| . (T4"5j) and (f4"oT) . we obtain the solution of the system with the 
initial data (xo,^o) 5 z o > 0, 



x(t) = t + iarccot M$ coth (f3(t)) 



t(t) = z exp /„' ^coth^( s) )^ ds 



i+ 



£0 coth(/J(s))) 



or 



ir(£) = i + iarccot &q tanh (/3(i)) 



*(t) = Z exp /„* Wanh(/3(s)) , 

V 1+ Ao tanh(/3(s))) y 



(49) 



(50) 



The graphs of the parametric curves in (|49j) , (|50j) are drawn in figure below 




Figure 6. Particle trajectory, in the irrotational case, for < c < 2 

As is shown in Figure 6, the coordinate z increases indefinitely in time, hence, 
the solutions (j49|) . (f50|) . are not physically acceptable. 
Thus, one gets the following theorem: 

Theorem 1. 

In the case that the underlying uniform current is moving in the same 
direction as an irrotational shallow water wave and the strength of the current 
is bigger than 2, then the particle trajectories beneath the wave are undulating 
paths to the right (see Figure 5). 

In the case that the underlying uniform current is moving in the opposite 
direction as an irrotational shallow water wave and the strength of the current 
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is smaller than -1, then the particle trajectories beneath the wave are undulating 
paths to the left (see Figure 3). If the strength of the adverse current is bigger 
than -1, then the particle trajectories are loops with forward drift (see Figure 4)- 
In the case of no underlying current and in the case that the underlying 
uniform current is moving in the same direction as an irrotational shallow water 
wave with the strength of the current smaller than 2, the particle trajectories 
obtained (see Figure 2 and Figure 6) are not physically acceptable. In these 
cases it seems necessary to study the full nonlinear problem. 

5.2 The case of a constant vorticity flow 

Making the Ansatz 

f(x-t)=cos(2ir(x-t)) (51) 
from (j3TJ)) . the differential system (|3"Tj) becomes 



§ = cos(2^(x - i)) + ^f^z + Co 



|| = 27rzsin(27r(x-t)) 



(52) 



Notice that the constant cq is the average of the horizontal fluid velocity on the 
bottom over any horizontal segment of length 1, that is, 

1 f x+1 

Co = — / u(s,0,t)ds, (53) 

1 .If 



The right-hand side of the differential system (|52|) is smooth and bounded, 
therefore, the unique solution of the Cauchy problem with initial data (xq, Zq) 
is defined globally in time. 

To study the exact solution of the system (|52"|) it is more convenient to re- 
write it in the following moving frame 

X = 2ir(x-t), Z = z (54) 

This transformation yields 



% = 2vrcos(X) + 2ir^f^-Z + 2tt(c - 1) 



H = 2nZsm{X) 



(55) 



Let us now investigate the differential system (1551) 

Differentiating with respect to t the first equation in (|55j) and taking into account 
(|55|) . this first equation becomes 

d 2 X 

— = 4tt 2 sin(X) [1 cq cospQ] (56) 
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Like in the irrotational case (see [12]), we use the following substitution (see 
[21], 1.76, page 308) 



cot 



X 



y , sin(X) = ^ V , cos(X) = ^ — \ , <1X = 



2 J ya + 1' ^ ' y 2 + l 

In the new variable, the equation (|56j) takes the form 



if + i 



oPy 2y / dy 



8^ 2 y 



47T 2 C £/ = 



dy 
(57) 

(58) 

(59) 
(60) 



with y(t) satisfying the ordinary differential equation (j58j) . 

From (|5l|) . ([5U1) and (|6U| . we obtain the solution of the system (|5"2"j) with the 

initial data (xq,zq): 



dt 2 y 2 + l\ dt J y 2 + l 
The solution of the system (j5"5j) has then the following expression: 

X(t) = 2arccot (y(t)), 



x(t) = t + iarccot (y(t)) 



z(t) = z exp\ 



(61) 



ds 



where y(t) satishes the ordinary differential equation (f5"5| with the following 
initial conditions: 



y(0) = cot(7ra ) 



dy 
dt 



(0)=7T 



9 



zo + co) (cot 2 (7ra ) 



(62) 



obtained from ((57)1 and the first equation in (|55l) . 

5.2.1 The case c = 

In this case the equation (|58|l writes as 



d 2 y 2y 



dt 2 y 2 + 1 



f 1 - 



For 



dy 
df 



7^ ±27r y 7^ ±27rf + a, with a constant , 



(63) 



(64) 



5 Particle trajectories beneath linearized shallow water waves 



15 



we put the equation (|63p into the following form 



-,y = 2 i » ( 65 ) 



(y)2_ 47r 2* y 2 + y 

where y' := % y" := gf. 

We integrate in (|65)l and we get 

(y') 2 - 4tt 2 = A\tf + If {y') 2 = A 2 (y 2 + l) 2 + 4vr 2 , (66) 

A being an integration constant, which, taking into account (1621) . can be ex- 
pressed function of the initial data (xq, zq) and of the constant vorticity ojq. 
The solution of the equation (|6"6l involves an elliptic integral of first kind, that 
is, 

±f- dv =t (67) 

The elliptic integral of first kind from (|67p may by reduced to Legendre's normal 
form. In order to do this we first consider the substitution 

y 2 = w (68) 

The left hand side in (l67l) becomes 



, dy _ , 1 f dw 



VA 2 (y 2 + l) 2 +4 7 r 2 2\A\ J ^ {u]2 ~ 2yj + " + 4£ } 
We introduce now the variable <p by (see 26 Ch. VI, §4, page 603) 



"'= \/ 1 + 5 tan2 2 (69) 



and we get 



, 2 4^ 2 \ / / ^r 2 \ , , 92 ,tan 2 f 



47T 2 tanf 
^ = V 1+ A 2 "^f^ 

where the constant < k 2 < 1 is given by 

f 




Therefore we obtain the Legendre normal form of the integral in (1671) , that is 

dip 

2\A\(1 + ^)J Vl-^sinV 
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Taking into account (jB"6")) . the derivatives of x(t) and z(t) from (|0 1 [) with 
respect to t, have the expressions 



, _ x(a 2 +i) T y^(^+i)^+4^ 

x W - 7r(y 2 +l) 

/ x (71) 



where y = y(i) is given implicitly by (|67|l . 

For the alternative with "+" in the expression (|71l) of x'(t), we obtain that 
x'(t) > 0, for all t. The sign of z'(i) depends on y(t). For j/(t) < we get 
z'(t) < 0, and for y(t) > we get z'(t) > 0. Concerning the zeros of the function 
y(t) we can tell more looking at the Legendre normal form of the integral in 
(j6"T|) . that is, at the expression (1701) . The inverse of this integral is, excepting 
a constant factor, the Jacobian elliptic function sn.(i):=sin.<£> (see, for example, 
Thus, taking into account (IU5|) . (|69p and the formula tan ^ = jj^^- , we 
get that the zeros of the function y(t) are the zeros of the periodic Jacobian 
elliptic function sn(t) (for the graph of sn(i) see, for example, [2], Figure 10, 
page 26). Therefore, the particle trajectory looks like in Figure 7 (a). 
Let us consider now the alternative with "-" in the expression ([71]) of x'(t). 
If 

\A\ > tt (72) 

then we have x'(t) < 0, for all t. The sign of z'{t) depends on y(t): for y(t) < 
we have z'(t) < 0, and for y(t) > we have z'(t) > 0. The function y(t) being 
periodic and having multiple zeros, the particle trajectory looks like in Figure 
7(c). 
If 

\A\ < tt (73) 
then the zeros of x'(t) are obtained for 



y(t)=±J-= J^ 1 (74) 

The sign of z'(t) depends on y(t) as above, that is, z'(t) < for y(t) < 0, and 
z'(t) > for y(t) > 0. Summing up, 

for y(t) € (-co, V J^_ A2 ~) we have x'(t) > 0, z'(<) < 0, 
for y(t) G {-^ v Ji 7T _ A2 ~ 0) we have x'(i) < 0, z'(t) < 0, 
for y(t) G (0, V J"L A 2 ~ 1) we have a/(i) < 0, z'(t) > 0, 
for G (^J v J*_ A i -A, oo) we have sc'(i) > 0, z'(t) > 0. 

If the periodic function y(t) has values in each interval from above, then the 
particle trajectory looks like in Figure 7 (b). 

Theorem 2. 
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As periodic waves propagate on the water's free surface of a constant vor- 
ticity shallow water flow over a fiat bed, with the average of the horizontal fluid 
velocity on the bottom over any horizontal segment of length 1 equals zero, there 
does not exist only a single pattern for all the water particles. The particle paths 
are not closed, and depending on the relation between the initial data (£n,Zn) 
and the constant vorticity ujq, some particle trajectories are undulating curves 
to the right (see Figure 7 (a)), or to the left (see Figure 7 (c)), others are loops 
with forward drift (see Figure 7 (b)). 



5.2.2 The case c ^ 

In this case we transform the equation ([58)1. by (see [21], 6.45, page 551) 

»M =-(§)'. < 75 » 

into the following equation 

2 dy y +1 y + 1 

The homogeneous equation: 

1 dw 2y 



2 dy ?/ 2 + l 

has the solution 



w{y) = (77) 



w h {y) = B(y 2 + l) 2 (78) 

where B is an integration constant. By the method of variation of constants, 
the general solution of the non-homogeneous equation (1761) is given by 

w(y)=B(y)(y 2 + l) 2 (79) 

where B(y) is a continuous function which satisfies the equation 

dB 87r 2 c y 16-7r 2 y 

rfy = (y 2 + i) 2 ~ (y 2 + 1) 3 ( } 

The solution of the equation (f80| is 

„, \ 47T 2 Co 47T 2 

where C is constant. Therefore, the solution of the non-homogeneous equation 
(j76|) has the expression 

w(y) = C(y 2 + l) 2 - 4^ 2 c (j/ 2 + 1) + 4tt 2 (82) 
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Taking into account (fT5j) , we now obtain, instead of the equation the 
following equation 



(y'Y = C(y 2 + I) 2 - A^cotf + 1) + 4^ 
where y' := and the following condition has to be satisfied 



(83) 



2 2 , n 

7T Cq < C 



(84) 



in order to have the right hand side in (|83[) bigger then zero for any y. 
The solution of the equation ((53)) involves an elliptic integral of first kind 



± 



/ 



dy 



^fC{y 2 + l) 2 - 47r 2 c (y 2 + 1) + 4tt 2 



(85) 



The elliptic integral of first kind from (|85l) may by reduced to Legendre's normal 
form. In order to do this we first consider the substitution 

y 2 = w 



The left hand side in (1851) becomes 



± 



±1 

2 



dy 



v/C(y 2 + l) 2 - 4ir 2 c {y 2 + 1) + 4tt 2 
dw 



Cw 



w 2 + w (2 



4ir 2 c 
C 



1 -I- 47T 2 47T 2 Cq 

1 ^ c c 



We introduce now the variable ip by (see [26] Ch. VI, §4, page 603) 



and we get 



a) 2- 



47T 2 C 



47T 2 47T 2 C 



4tt 2 47r 2 c \ , 
C C / 



fc 2 sin 2 y?) 



tan 2 f 
co?f" 



, 4tt 2 4ir 2 c tanf 
rfW=A/l + — - C cos 2 * ^ 



where the constant < k 2 < 1 is given by 



; 1- 



C - 2n 2 c 



CJ1 + ^-^ 
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Therefore we obtain the Legendre normal form of the integral in (|85[). that is, 



2^(1 + ^-^0)^ Vl - * 2 sin 2 <p 



Taking into account ([83]). the derivatives of a;(i) and z(t) from ([61]) with 
respect to t, have the expressions 



. _ 7r(y 2 + l) Ty /C(y 2 + l) 2 -4:^c (y^ + l)+47r^ 
X \ l > ~ 7r(a 2 + l) 



(87) 



z 



where y = y(t) is given implicitly by 
For the alternative with "+" in the expression (|87l) of x'(t), we obtain that 
x'{t) > 0, for all t. The sign of z'it) depends on yit). For y(t) < we get 
z'(t) < 0, and for y(t) > we get z'it) > 0. Taking into account the Legendre 
normal form of the integral in (|85[) . that is, the expression in (1861) . and using 
the same arguments as in the case en = on the page 14, we get that the zeros 
of y(t) are the zeros of the periodic Jacobian elliptic function sn(t). Therefore, 
in this case the particle trajectory looks like in Figure 7 (a). 
Let us consider now the alternative with "-" in the expression ([87)) of x'(t). The 
zeros of x'(t) are obtained by solving the equation 

(tt 2 - C)(y 2 + l) 2 + 47r 2 c (y 2 + 1) - 4tt 2 = (88) 

The discriminant of the quadratic equation in W 

(tt 2 - C)W 2 + 4tt 2 c W -At: 2 = (89) 



is 

A = 16tt 2 (t7 2 Cq + tt 2 — C) (90) 

I) If 

C>n 2 c 2 +n 2 (91) 

then A < and tt 2 - C < 0. This yields x'(t) < 0, for all t. The sign of z'(t) 
depends on y(t) as was mentioned above. The function y(t) being periodic and 
having multiple zeros, the particle trajectory looks like in Figure 7 (c). 

II) If 

C < TT 2 C 2 + TT 2 (92) 



then A > and the equation (|89|) has two real solutions. 
II a) If in addition to the condition (|92|) we have 



tt 2 - C > 



(93) 
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then, one of the solution of the equation ([55)1 is positive and the other one is 
negative. We denote the positive solution of the equation (|89| by 



TX Z — C 

Further, if W\ — 1 < 0, then the sign of x'(t) is given by it 2 — C, that is, in this 
case x'(t) > 0, for all t. Taking into account that z'(t) > for y(t) > 0, and 
z'(t) < for y(t) < 0, the particle trajectory looks now like in 7(a). 
If Wi — 1 > 0, then the zeros of x'(t) are obtained for 



y(t)=± y /W 1 -l (95) 

The sign of z'(t) depends on y(t) as above. Thus, 
for y(t) e (-oo, -y/W[ - 1) we have x'(t) > 0, z'(t) < 0, 
for y(t) e (-V Wi - 1, 0) we have x'(t) < 0, z'(t) < 0, 
for y(t) e (0, VWi - 1) we have x'(t) < 0, z'{t) > 0, 
for e (VWi - 1, oo) we have x'{t) > 0, z'(t) > 0. 

If the periodic function y(t) has values in each interval from above, then the 
particle trajectory looks like in Figure 7 (b). 

II b) If in addition to the condition (|92[) we have 

TT 2 - C < (96) 



then, the solutions of the equation (|89|) are both negative or both positive. 
If 

c < (97) 

then both solutions of (|89| are negative. Therefore, the sign of x'(t) is given by 
7r 2 — C, that is, in this case x'(t) < 0, for all t. How z'(t) < for y(t) < and 
z'(t) > for y(t) > 0, the particle trajectory looks like in 7(c). 
If 

c > (98) 
then both solutions of (|59"j) are positive. We denote these positive solutions by 



Wi 



-2tt 2 c + 27r v / 7 r 2 c§ + tt 2 - C 



W 2 := ~ 27r2c " - ± gEg (99) 

We observe that W\ <W 2 - 

Further, if W\ — 1 < and W 2 — 1 < then the sign of x'(t) is given by n 2 — C, 
that is, in this case x'(t) < 0, for all t. Therefore, the particle trajectory looks 
like in 7(c). 

If W\ — 1 < and W 2 — 1 > then the zeros of x'(t) are obtained for 



y(t) = ±^W 2 - 1 



(100) 
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The sign of z'(t) depends on y(t) as above. Thus, taking into account 
for y(t) G (-o o,-VW 2 - 1) we get x'(t) < 0, z'{t) < 0, 
for y(t) G (-V W 2 - 1, 0) we get £c'(t) > 0, < 0, 
for y(t) € (0, -v/Wg - 1) we get a/(t) > 0, z'(*) > 0, 
for t/(t) G - 1, 00) we get x'{t) < 0, z'(i) > 0. 

If the periodic function y(t) has values in each interval from above, then the 
particle trajectory looks like in Figure 7 (d). 

If W\ — 1 > and W2 — 1 > then the zeros of x' (t) are obtained for 

y(t) = ±y/Wi - 1, y = ±y/W 2 -l (101) 

The sign of z'(t) depends on y(t) as above. Thus, taking into account (|M|) . 

for y(t) G (-0 0, —y/ W 2 - 1) we ha ve x'(t) < 0, z'{t) < 0, 

for y(t) e (- VW 2 - 1 , -y/Wi - 1) we have sc'(t) > 0, z'(t) < 0, 

for G (--y/ Wi - 1, 0) we have x'(t) < 0, z'(t) < 0, 

for y(t) € (0, -y/Wi - 1 ) we hav e x'(t) < 0, z'(t) > 0, 

for G (V lVi - 1 , VW 2 - 1) we have cc'(t) > 0, z'(t) > 0, 

for y(t) G (VW2 - 1, 00) we have x'(t) < 0, > 0. 

If the periodic function y(t) has values in each interval from above, then the 
particle trajectory has the peculiar shape in Figure 7 (e). 

Theorem 3. 

As periodic waves propagate on the water's free surface of a constant vorticity 
shallow water flow over aflat bed, with the average of the horizontal fluid velocity 
on the bottom over any horizontal segment of length 1 different from zero, there 
does not exist only a single pattern for all the water particles. The particle paths 
are not closed, and depending on the relation between the initial data (xq,Zq) 
and the constant vorticity ujq, some particle trajectories are undulating curves 
to the right (see Figure 7 (a)), or to the left (see Figure 7 (c)), others are loops 
with forward drift (see Figure 7 (b)), or with backward drift (See Figure 7 (d)), 
others can follow peculiar shapes (see Figure 7 (e)). 
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(a) 



(b) 




(c) 



(d) 




(e) 




Figure 7. Particle trajectory in the constant vorticity case 
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